Ubiquity of optical activity in planar metamaterial scatterers 
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Recently it was discovered that periodic lattices of metamaterial scatterers show optical activity, even if 
the scatterers or lattice show no 2D or 3D chirality, if the illumination breaks symmetry. In this Letter we 
demonstrate that such 'pseudo-chirality' is intrinsic to any single planar metamaterial scatterer and in fact has 
a well-defined value at a universal bound. We argue that in any circuit model, a nonzero electric and magnetic 
polarizability derived from a single resonance automatically imply strong bianisotropy, i.e., magneto-electric 
cross polarizability at the universal bound set by energy conservation. We confirm our claim by extracting 
polarizability tensors and cross sections for handed excitation from transmission measurements on near-infrared 
split ring arrays, and electrodynamic simulations for diverse metamaterial scatterers. 



Many historical debates on how to describe the effec- 
tive electrodynamic response of media composed of sub- 
wavelength building blocks currently acquire new relevance 
in nano-optics. Initiated by the works of Veselago and 
Pendry U efforts are focused on manipulating effective 
medium parameters in nanostructured media. On the one 
hand, the drive for arbitrary e and /i is generated by the idea 
that light fields can be arbitrarily reshaped by conformal trans- 
formations, provided we can create arbitrary constitutive ten- 
sors Isj-IH]. On the other hand, a convergence with plasmon- 
ics has led to the realization that subwavelength scatterers 
mimic and even greatly enhance rich scattering phenomena 
known from molecular matter. For example, resonantly in- 
duced optical magnetism in 2D and 3D chiral metal nano- 
objects results in giant circular birefringence, optical rotatory 
power, broadband optical activity, and circular dichroism in 
frequency ranges from microwave, mid-IR, near IR to even 
visible frequencies ll6[- fl4ll . The fact that strong optical ac- 
tivity is easily attained using chiral subwavelength scatterers 
is promising for many applications such as broadband op- 
tical components, as well as providing excellent candidates 
for achieving negative refraction jisll . or repulsive Casimir 
forces lfl6ll . Moreover, the promise of enhancing detection 
of molecular chirality via enhanced chirality in the excitation 
field, is expected to be of large importance for, e. g., d iscrimi- 
nation of enantiomers in biology or medicine il7l420ll . 



A question of essential importance is how to control the 
optical activity of a single building block, i.e., have indepen- 
dent control over the degree of magnetic response, electric re- 
sponse and magneto-electric cross coupling or 'bianisotropy' 
whereby incident electric (magnetic) fields cause a magnetic 
(electric) material polarization in a single building block 12111 . 
For instance, in attempts to reach negative indices, researchers 
soon found that the archetypical split ring resonator (SRR) 
has a magneto-electric response that is undesirable, yet diffi- 
cult to remove without also losing the magnetic response 12211 . 
Completely opposite to the desire to remove this bianisotropy, 
it has also been realized that all applications exploiting op- 
tical activity benefit from strong magneto-electric coupling. 
Currently it is unclear if there exists any universal bound to 
which optical activity can be benchmarked, or conversely, if it 
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FIG. 1 . a) Any scatterer a with nonzero electric and magnetic polar- 
izability shows oblique incidence optical activity, with transparency 
for one handedness of incident light at off-angles, and maximum ex- 
tinction when the incident beam is rotated by 90°. At normal in- 
cidence, the scatterer shows no optical activity, b) Common planar 
scatterers for which we verify optical activity and bianisotropy: (1) 
scanning electron micrograph of 230x30x30 nm Au SRRs. Struc- 
tures (2)-(6): n particles of varying arm length. Structure (7) model 
for SRR in (1). Structure (8,9,10): double split ring and double gap 
ring m. 



is at all possible to avoid bianisotropy without also losing the 
magnetic response 12311 . In this Letter, we discuss precisely 
such a universal bound for magneto-electric coupling for sin- 
gle scatterers, disentangled from any lattice properties. We 
claim that Onsager's relations constrain optical activity to al- 
ways be at this maximum bound for any dipole scatterer based 
on planar circuit designs, independent of geometrical chiral- 
ity. Our claim is supported by measurements on SRRs at tele- 
com wavelengths and rigorous fuU wave calculations ll24ll in 
which we retrieve cross sections and polarizabilities for vari- 
ous metamaterial scatterers (see Fig.[na,b)). 



The central quantity in this Letter is the polarizability tensor 
that quantifies the magnetic response, electric response and 
magneto-electric cross coupling (bianisotropy) intrinsic to a 
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FIG. 2. Transmission spectra from a periodic square array of 
230x230 nm SRR with ci=530 nm. Tiie spectra were taken as a 
function of angle of incidence, wiiere daslied curves denote negative 
angles, and solid curves positive angles with respect to the sample 
normal. (a,c) and (b,d) are transmission spectra shown for right- and 
left-handed circularly polarized illumination. Inset in (a) resp. (c) 
shows the the sample rotation axis for (a,b) resp. (c,d), taking the 
incident fc-vector as pointing through the paper 



single metamaterial building block according to 112 IL 12311 
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For molecules, optical activity is due to weak cross coupling, 
i.e., a perturbative acH w 10~^ asE, while an ~ 0. In 
contrast, the paradigm of metamaterials is that a single scat- 
terer acquires a magnetic dipole moment m at least compa- 
rable to the electric moment p, with ue, an, and possibly 
ac of the same order, which all derive from a single res- 
onance II25II . In order to quantify the polarizability for the 
canonical SRR, we performed transmission measurements as 
well as full-wave calculations. For the experiments we fabri- 
cated Au SRRs arranged in square arrays on glass substrates 
by electron beam lithography (e-beam), resonant at telecom 
wavelengths |26, 2^. Fig.fTJb) shows a scanning electron mi- 
crograph (SEM) of a SRR array with 530 nm lattice spa cing , 
which is so dilute that coupling between SRRs is small 12711 . 
yet so dense that no grating diffraction occurs. Each SRR 
measures 230 x 230 x 30 nm, with a gap between the arms that 
is 100 nm wide and 145 nm deep. We record transmission by 
illuminating the sample with a narrow band of frequencies at 
a time, selected from a supercontinuum laser (Fianium), using 
an acousto-optical tunable filter (Crystal Technologies) with a 
bandwidth of 1-2 nm 1281]. The beam is chopped for lock-in 
detection on an InGaAs photodiode. We polarize the incident 
beam using a broadband quarter-wave plate, to provide cir- 
cularly polarized excitation. We weakly focus the beam onto 
the sample (/=100 mm). Light is collected with a low NA 



collection lens (/=20 mm), and passed through a telescope 
and pinhole to ensure spatial selection from within a 200x200 
finP e-beam write field, as monitored by an InGaAs camera. 
A motorized rotation stage allows transmission measurements 
versus incident angle relative to the sample normal. 

Fig. |2] shows transmission versus wavelength for left and 
right handed circularly polarized incident light, for incidence 
angles from -50° to H-50°. Fig. |2](a) shows data when the 
angle is varied from normal incidence by rotating the SRRs 
around their mirror axis y. At normal incidence, the magnetic 
LC resonance is evident around 1600 nm wavelength as a min- 
imum in transmission. As opposed to the deep minima usually 
reported for linear, x-polarized transmission {E along the gap) 
of dense arrays, the transmission dip is shallow since our lat- 
tice is dilute and the LC resonance is associated only with E^ 
and Hz, and completely transparent for Ey. As the incidence 
angle is moved away from the normal, the excitation also of- 
fers Hz as a driving field, a quarter wave out of phase with 
Ex- A very clear asymmetry around the normal develops. For 
right-handed light the transmission minimum becomes con- 
tinuously shallower towards negative angles, and the sample 
is nearly transparent for —50°. In contrast, the transmission 
minimum significantly deepens from 28% to 75% when going 
to large positive angles. The asymmetric behavior with inci- 
dence angle is mirrored for opposite handedness (Fig. |2jb), 
consistent with oblique incidence optical activity. For linear 
polarization the transmission is symmetric around normal in- 
cidence (not shown). 

The fact that optical activity is symmetry-allowed even for 
lattices containing 2D non-chiral objects aligned with the lat- 
tice symmetry, was already reported by Plum et al. (2^], who 
coined this 'extrinsic 3D chirality'. In contrast to symmetry 
arguments that only distinguish between allowed and forbid- 
den effects without quantifying the strength of optical activity, 
it is the express aim in this Letter to ascertain what the single 
element polarizability is that leads to the strong optical ac- 
tivity. We exclude the array structure factor as the cause of 
handed behavior lisoll , as the optical activity disappears when 
we rotate the SRRs by 90° in the sample plane (Fig.|2](c) and 
(d)). We hence conclude that the single SRR polarizability 
must contain the strong 'pseudo-chiralilty' that is expressed 
as huge circular dichroism contrast in the extinction cross sec- 
tion, despite SRRs being neither 2D nor 3D chiral. Qualita- 
tively, the LC description of a single SRR indeed contains 
optical activity under oblique incidence. Charge motion is set 
by g = {iujL + R+l/C)~^[iuJiJoAHz+E.j:/t],where Listhe 
inductance, C the capacitance, R the Ohmic resistance, t the 
capacitor plate gap and A the enclosed area. Full transparency 
despite the presence of suitable driving E^ along the gap and 
Hz through the split ring occurs when iujuqAHz = —E^jt. 
Conversely, optimum driving of a SRR benefits from an oppo- 
site quarter wave phase difference between E.j^ and H^ so that 
^uj^uqAHz + i^x/t] has maximum magnitude. Circular po- 
larization at oblique incidence provides the required quarter 
wave phase difference between E^ and H^. 

We quantify the polarizability tensor from the data by ana- 
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FIG. 3. Circles: effective extinction per SRR from transmission 
data. Solid line: single scatterer extinction cross section expected 
in a dipole model. Dashed line: lattice sum calculation for a square 
array with pitch d=530 nm of magneto-electric dipoles. 



lyzing the effective extinction cross section per SRR defined 
as (T = {1 — T)(P, where d is the lattice pitch and T is the min- 
imum in transmission 12711 . Fig.[3]shows that this effective ex- 
tinction cross section varies between 0.07 and 0.16 /im^ as the 
angle is swept from ±50° to ^50° (mirrored dependence for 
opposite handedness). For a single magneto-electric dipole 
scatterer Ref. lEsll predicts that the extinction cross section 
generally depends on angle 9 as 



trA^) =<^- + - + cos{2{e ± 9o))]/2. (2) 

Measurements on a single object would provide the elec- 
trodynamic Isill as through the normal incidence extinc- 
tion ct_r,l(0) = 2TTklmaE, while the maximum and min- 
imum attained extinction a± encode electrodynamic po- 
larizability eigenvalues via a± = TTkIm{aE + an ± 
\J {oiE — cen)^ + 4q;^). Such a fit of the single object ex- 
tinction to the measured effective extinction would provide 
aE = 4:.1V, an = i-QV and ac ~ lAV expressed in units 
of the geometrical volume of the SRR (V = 0.0012 ^m^). 
However, in a lattice of SRRs, the response is modified by co- 
herences such that (p, m)-^ = 1/[q:~-'^ — fJ](E,H)-^, where 
a lattice sum Green function Q renormalizes the polarizabil- 
ity 11321]. We calculate lattice transmission by rigorous electro- 
dynamic lattice sums involving all multiple- scattering interac- 
tions between SRRs t32<l . Consistent with our data, the calcu- 
lated transmission shows strong optical activity under oblique 
incidence. We extract as = 6AV, an = 0.9y , ac = 2.1V 
at A=1600 nm from a comparison to data, highlighting that the 
response of SRR arrays is consistent with remarkably strong 
maximum magneto-electric cross coupling. 

In Ref. 123'] we analyzed how electrodynamic scatterers 
with arbitrary polarizabilities of the form in Eq. ([T]i scat- 
ter. In that work, we realized that once one applies the op- 
tical theorem to a planar scatterer (in-plane p, out-of-plane 



m), ac < y/oEOH appears as the maximum value that 
ac - the crosscoupling after taking a common resonant fre- 
quency factor out of Eq. (HJ \2% - can possibly attain to 
avoid violation of energy conservation. Here we claim that 
any planar circuit-derived scatterer is necessarily exactly at 
this upper bound, i.e., at maximum cross coupling. To prove 
this assertion we analyze a generic model for the polariz- 
ability of a planar scatterer under two general assumptions: 
(1) a linear response and (2) that an electric and magnetic 
dipole response originate from the same equation of motion 
for charge q moving through the scatterer. Linear response 
impHes q = Ce{^^)E + Ch{'-o)H, where E (H) is in the 
plane (perpendicular to the plane) of the scatterer Since p 
and m both derive from the same charge motion, p = Apq and 
m = A,nq = iLj-!Am{i^)q, where Ap and Am are geometry- 
dependent constants. One now finds the electrostatic circuit 
polarizability as 



ao 



ApCe{uj) 

iujArnCE{(^) 



ApCniuj) 
iujA„iCHiuj) 



(3) 



For reciprocal materials, Onsager's relations constrain aE 
and aH to be symmetric, as well as requiring ^pCi/(w) = 
—iujAmCE{}jj)- Taking out a common frequency factor 
C{ijj) oc Ce{uj) that describes the circuit resonance, one finds 
that ao always take the form ll25[l 
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-lujy/aEaH 
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(4) 



The surprise is that Onsager constraints leave no freedom to 
choose the off-diagonal coupling ac- Any planar circuit ele- 
ment is cross coupled, with cross coupling ac = \/aF,aH- 
Combining this finding with our result from Ref. ll23ll we 
conclude that any planar circuit-derived scatterer is not just 
cross coupled, but that this coupling is at the maximum cross 
coupling limit. Maximum cross coupling means one vanish- 
ing eigenpolarizability a- = 0, hence complete transparency 
of the scatterer for one handedness under oblique incidence, 
which means huge optical activity contrast. 

Based on our experiment, we can now assess whether the 
strong cross coupling in real scatterers is indeed close to the 
predicted maximum. From the polarizability we extracted 
from the very strong circular polarization contrast in extinc- 
tion observed for split rings in Figure[3]we indeed find almost 
maximum cross coupling, since ac ~ O.SS^a^a//. Fur- 
thermore, we use full-wave simulations to examine the po- 
larizabiUty, and pseudo-chirality in extinction of many scat- 
terers. We use 3D Surface Integral Equation (SIE) calcu- 
lations 12411 . to obtain full-wave solutions for archetypical 
metamaterial scatterers including SRRs, Omega particles with 
straight legs of different length, double SRRs and double-gap 
rings as shown in Fig.lHb). We calculate scattering cross sec- 
tions and polarizability tensors independently from each other 
To extract the polarizability, we excite the same scatterer with 
six linearly independent illumination conditions, obtained as 
counter-propagating linearly polarized beams set in (out of) 
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FIG. 4. Master diagrams summarizing optical activity and bi- 
anisotropy mapped as a function of ^ = {as — aH)/{{aE + an) 
and 77 = [otc) / {{ctE + olh) ■ All structures we tested (data-points, 
numbered as in Fig. [Tfb)) are close to the locus of maximum cross 
coupling (ellipse), except (8). The color scale shows optical activity 
contrast in the dipole approximation (color scale) and for tested 
structures (dots). Panel (b) is a 3D representation of (a). 



phase to yield just electric (magnetic) Cartesian excitation. 
We project the calculated scattered E-field evaluated on a 
spherical surface around the scatterer on vector spherical har- 
monics to retrieve p and m l33ll . As a consistency check on 
the polarizability retrieved by matrix inversion we verify that 
the Onsager constraints are satisfied, which are not a priori as- 
sumptions in the retrieval. We summarize results for all scat- 
terers in a 'master plot' that allows comparison independent of 
scatterer size. The scatterers are shown in Fig.[T](b). As a first 
dimensionless variable we use f = [as — aH)/{aE + 01 h), 
which equals ±1 for purely electric (magnetic) scatterers, and 
for equal electric and magnetic polarizability. As a dimen- 
sionless second variable we take the normalized cross cou- 
pling 1] — ac /{aE + c(h)- The locus of maximum cross cou- 
pling is the ellipse r; 



^/l - Fig. m shows that most 

metamaterial scatterers we analyzed have ^ well away from 1, 
indicating significant magnetic polarizability. Furthermore all 
particles are essentially on the locus of maximum cross cou- 
pling, confirming our claim that bianisotropy is ubiquitous. 

As third axis for the master plot we use a measure for op- 
tical activity in scattering. All scatterers we simulated show 
an angular dependence of the scattering cross section of the 
form in Eq. (|2]i. The dimensionless parameter '5 = \aR — 
ctl|/ (cTfl + cTL) evaluated at 45° incidence angle quantifies the 
maximum attained difference in extinction {an — (Jl| (maxi- 
mal always at 45°) normalized to (twice) the angle-averaged 
extinction cross section (7+ + (7_. Fig. |4] shows versus ^ 
and 77 as predicted by point scattering theory. Evidently, opti- 
cal activity is expected to be absent for zero cross coupling, 
and to increase monotonically as cross coupling increases. 
Very strong contrast in extinction per-building block is ex- 
pected along most of the locus of maximum cross coupling. 



vanishing only for purely electric, and purely magnetic dipole 
scatterers (^ = ±1). The full-wave simulations show that all 
the commonly used metamaterial scatterers exhibit strong op- 
tical activity in surprisingly good agreement with the dipole 
model given that the circuit approximation, and the neglect of 
multipoles and retardation in Eq. |4] are very coarse assump- 
tions. Freedom to deviate significantly from the dipole model 
requires multiple overlapping resonances in a single scatterer. 
Indeed, the most noted deviations occur for the object (8,9) 
which has two hybridized resonances of separate parts. Ear- 
lier findings based on symmetry arguments proposed that ex- 
trinsic 3D chirality requires loss 12911 . We find that optical 
activity is in fact ubiquitous for planar magneto-electric scat- 
terers, irrespective of absorption. The cancellation of optical 
activity for zero absorption noted by does not occur in 
a but occurs in special cases where observables are subject 
to additional symmetries, such as wave vector conservation in 
non-diffracting periodic systems. 

To conclude, we have shown that planar metamaterial scat- 
terers that rely on a single resonance to generate a simul- 
taneous electric and magnetic response are maximally bian- 
isotropic and strongly optically active, whether they exhibit 
geometrical chirality or not. Our findings have important im- 
plications for controlling bianisotropy independently of e and 
/i in metamaterials, since they imply that it is fundamentally 
impossible to independently control bianisotropy for single 
resonant objects. The only route to avoid bianisotropy in lat- 
tices of resonators is to use heterogeneous lattices that contain 
distinct, or multi-resonant elements (e.g., double-split rings 
in Fig. lU to independently generate e and /i, or to use lat- 
tices of effectively larger 'super-cells' with rotated copies of 
the same building block to cancel off-diagonal coupling. Our 
results also hold important promise for enhancing far-field or 
near-field chirality 11911 in scattering applications where it is 
desired. In general, since maximum cross coupling is ubiqui- 
tous, optical activity is a very robust phenomenon that is easily 
extended to, e.g., finite clusters, random assemblies, or multi- 
element antennas. For instance, we predict that one can create 
chiral variants of the plasmon Yagi-Uda antenna to generate 
or selectively enhance circularly polarized single emitters. 
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